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ABSTRACT 

Aims. Images obtained with the ESO VLT and FORS1 in [O in] 5007 on- and off-band, as well as r_Gunn filters, are analyzed to 
search for planetary nebula (PN) candidates. 

Methods. In the continuum-subtracted [O in] 5007 on-band images, a large number of emission-line regions were detected. We 
describe the criteria employed for distinguishing PN candidates from compact HII regions. 

Results. The most unambiguous discriminators for the two classes of nebulae are the sizes and properties of the ionizing stars. Based 
upon these criteria, we have found 20 PN candidates for which we present coordinates, nebular [O m] fluxes, and stellar magnitudes. 
The cumulative luminosity function for these PNe is discussed. A catalog of HII regions listing coordinates, nebular [O in] fluxes, 
stellar magnitudes, and other characteristics is also presented. We find that HII regions are rather concentrated towards the disk of the 
galaxy, while PNe are distributed also above and below this structure, consistent with their belonging to an older stellar population. 
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1. Introduction 

The dwarf galaxies of the Local Group (LG) and its vicin- 
ity allow detailed studies of their different component popula- 
tions, such as stars, planetary nebulae (PNe), and HII regions. 
Therefore, they provide strong constraints for stellar population 
studies and help in understanding how the elemental abundances 
of galaxies have changed over time. Dwarf irregulars, in partic- 
ular, allow a comparative analysis of the chemical properties of 
their stellar, PN, and HII region populations, which provides a 
probe of the chemical enrichment and star formation history of 
the galaxy from intermediate ages to the present day. In addition, 
through spectroscopic analysis of PNe and HII regions in differ- 
ent zones of these galaxies, it is possible to address the ques- 
tion of the large-scale homogeneity in the distribution of heavy 
elements. To proceed with these studies, the first step is to de- 
tect PNe and HII regions, distinguishing them from each other, 
which is a non-trivial matter in star-forming galaxies, before ob- 
taining spectrophotometry of the candidates. Spectroscopic anal- 
yses have already been done for several galaxies of the Local 
Group. For example, Stasiriska et al. (1998) present a compar- 
ative analysis of the PN populations in the Magellanic Clouds, 
M32, and the bulges of the Milky Way and M31. M33 was an- 
alyzed by Magrini et al. (2003), while Sextans A and B were 
studied by Magrini et al. (2005) and Kniazev et al. (2006). 
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Here and in a forthcoming paper, we propose studying the 
PN and HII region populations in NGC 3109 (DDO 236), a 
Magellanic-type irregular classified as a late-type spiral SB(s)m 
(de Vaucouleurs et al. 1991). NGC 3109 is the dominant mem- 
ber of the small group of metal-poor dwarf galaxies that also in- 
cludes Sextans A, Sextans B, and Antlia, just at the edge of the 
LG. Due to its distance, NGC 3109 is particularly well-suited to 
chemical evolution studies because it does not seem to have been 
affected by tidal effects from the massive galaxies of the LG, so 
it should not have lost a significant amount of gas. 

With an absolute magnitude Mb ~ -15.2 (see §4) and a 
mean oxygen abundance log O/H+12 ~ 7.75 (as derived from 
HII regions by Lee et al. 2003 and, recently, from blue super- 
giants by Evans et al. 2006), NGC 3109 is similar to the SMC in 
luminosity and chemical composition. The galaxy is seen almost 
edge-on, having an elongated shape aligned almost E-W, with an 
extension of about 17.4x3.5 arcmin. Studies of the stellar popu- 
lations in this galaxy have been made by, among others, Demers 
et al. (1985), who pointed out some traces of spiral structure. 
Sandage & Carlson (1988) reported light curves and periods for 
29 cepheids and estimated a distance modulus of 26.0 mag (1.59 
Mpc). The cepheids were revisited by Capaccioli et al. (1992), 
who derived a distance modulus of 25.5+0.16 (1.26+0.10 Mpc), 
and by Musella et al. (1997), who find a distance modulus of 
25.67+0.16 (1.36+0.10 Mpc). Based upon V and I band pho- 
tometry, Minniti et al. (1999) established the presence of a halo 
consisting of Population II stars and derived a distance modu- 
lus of 25.62+0.1 mag (1.33+0.10 Mpc) from a deep luminosity 
function of these stars. We adopt this value for the distance to 
NGC 3109 in the following sections. 
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NGC 3109 presents all of the usual indicators of robust star 
formation activity: numerous extended and compact HII regions, 
large ionized shells, super-shells, and filaments of low surface 
brightness. A catalog of extended HII regions was published by 
Hodge in 1969. The OB associations and the distribution of star- 
forming regions were analyzed by Bresolin et al. (1993), who 
also produced a catalog of HII regions. Hunter et al. (1993) ana- 
lyzed the shells and super-shells in this galaxy. 

Richer & McCall (1992) reported the discovery of seven PN 
candidates and ten HII regions in the central zone of NGC 3109 
and, from the brightest-PNe method, they derived a distance 
modulus of 25.96 (+0.1/-0.4) corresponding to 1.6 Mpc. From 
good-quality imaging Pena et al. (2006; and this paper) have 
found that most of the PN candidates found by Richer & McCall 
are more likely compact HII regions that are somewhat brighter 
than PNe, which explains why their distance modulus is overes- 
timated compared to the stellar studies. Leisy et al. (2005) report 
the discovery of 18 PN candidates in NGC 3109 and, from spec- 
troscopy, they find that PNe display a wide range of metallicities. 

Here, we present a deep survey of PNe and HII regions in 
NGC 3109 based on [O m] 5007 on- and off -band images, as 
well as rjGunn imaging covering almost the entire galaxy. The 
data were acquired with the ESO Very Large Telescope (VLT), 
equipped with the Focal reducer spectrograph FORS1. These 
images constitute the 'pre-imaging' observations for the spectro- 
scopic program ID 076.B-0166. They were obtained to identify 
candidates for a subsequent multi-object spectroscopy. From the 
continuum-subtracted [O in] 5007 images, we identified more 
than 40 HII regions and 20 PN candidates that are reported in 
this work. 

The VLT-FORS 1 multi-object spectroscopic run associated 
with the pre-imaging observations was performed on January 
27, 2006. Spectra of several HII regions and PN candidates, dis- 
tributed in different zones of the galaxy, were acquired. These 
data allowed us to confirm the PN nature of many PN candidates 
and will form the basis of a forthcoming paper (Pena et al., in 
preparation). Nonetheless the calibrated [O ra] 5007 fluxes are 
presented here, as they allow us to calibrate the [O m] on-band 
imaging and help to strengthen our arguments. In our next paper, 
we shall study the chemical abundances of nebulae and star for- 
mation history of the galaxy, along with the large-scale spatial 
distribution of the heavy elements. 

This paper is organized as follows: in §2 we present the 'pre- 
imaging' observations and discuss the criteria we use to dis- 
tinguish between planetary nebulae and compact HII regions. 
In §3, we present the properties of the detected PNe and HII 
regions, including positions, diameters, [O m] 5007 calibrated 
fluxes, nebular magnitudes, and V magnitudes for the central 
stars. The cumulative luminosity function for the PNe is calcu- 
lated and discussed in §4, and the conclusions presented in §5. 

2. Imaging: distinguishing between PNe and 
compact HII regions 

The wide field coverage of FORS 1 allowed us to observe almost 
the entire galaxy with only two frames. The observing program 
consisted of imaging with the filters [O m] 5007 on-band (filter 
name FILT_500_5+85, ^ e =5000 A, FWHM=50 A), [O m] off- 
band (filter name OIII/6000+52, Ac= 5105 A, FWHM=61 A) 
and rJGunn (hereafter r). The images were acquired in service 
mode on November 29, and December 1, 2005. 

Two fields were observed, centered at 10:03:19.8 -26:09:32 
(east field; ob ID 204405) and 10:02:54.5 -26:09:22 (west field; 



ob ID 204406). Each field has dimensions of 6.8x6.8 arcmin 
on the sky (they overlapped by about 2 arcmin in the E-W 
direction), with a plate scale of 0.2"/pix. For each field, four 
dithered images were acquired with each filter, each with ex- 
posure times of 94 s (on-band), 90 s (off -band), and 10 s (r). The 
images were reduced and calibrated with the normal procedures 
of the ESO pipeline. That is, reduced images are bias-subtracted, 
flat-fielded, and astrometrically calibrated. The astrometric pre- 
cision is 0.2 arcsec. 

Each set of four images was combined using irah3 rou- 
tines, allowing for cosmic ray removal. The resulting images are 
equivalent to a 6.2 min exposure for the [O m] on-band filter, 6 
min for the off -band filter, and 40 s for the r filter (these images 
will be available at the CDS). The on- and off-band combined 
images were subtracted (after background subtraction and nor- 
malization) to easily detect [O m] 5007 emission objects. 

No standard stars were observed for the photometric calibra- 
tion of the images. Nevertheless in §3 we show that a rough cal- 
ibration of the stellar magnitudes can be performed using some 
known stars from the field, and the [O m] on-band images can 
be calibrated using our spectroscopy. 

The continuum-subtracted [O m] 5007 images show numer- 
ous [O in] 5007-emitting objects. The 7 PN candidates and the 
10 HII regions reported by Richer & McCall (1992) were re- 
covered. At least four of the large shells studied by Hunter et 
al. (1993) are clearly visible; the others are not, probably be- 
cause they are low excitation nebulae that do not emit in [O ra] . 
Several of the HII regions reported by Hodge (1969) were also 
identified. In addition, a number of small diffuse HII regions ex- 
ist, as well as many unresolved emitting objects. Some of them 
coincide with the Ha emitting objects in the list from Bresolin 
et al. (1993), but the coordinates provided by these authors have 
uncertainties larger than 1", and sometimes it is difficult to iden- 
tify their objects with ours. Our objects are presented in Tables 
1 and 2, where we have included a column named "other ID" 
listing the names given by other authors for the same object. 

One of the most important issues in studying extragalactic 
PNe is to distinguish true PNe from compact HII regions and, 
of course, to avoid contamination with background objects. To 
identify PNe in our imaging, we used the criteria described be- 
low. Criteria (a) and (b) are very similar to those employed by 
Soffner et al. (1996) in their work on NGC 300. 

a. Object size 

At the distance of NGC 3109 (about 1.33 Mpc), one arcsec 
corresponds to 6.4 pc. Given our pixel size of 0.2 "/pix and our 
FWHM of about 3.2 pix for stellar objects, PNe (which typically 
have diameters smaller than 1 pc) should appear as point-like ob- 
jects. Thus, we selected the unresolved [O m]-emitting objects. 
These nebulae can be either PN candidates or unresolved HII 
regions. 

b. The central star 

HII regions and PNe with significant [O m] 5007 emission 
are ionized by stars with T e fr > 35,000 K (e.g., Stasihska 1990). 
In the case of HII regions, this implies a central star of spectral 
type 08V or earlier (or 07 if the luminosity class is III or la). 
Typically, such stars have absolute visual magnitudes of My ~ 
-4.3 to -4.6 or brighter (Martins et al. 2005; Vacca et al. 1996). 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract to the National Science 
Foundation. 
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On the other hand, the central stars of galactic PNe with T e ff > 
35,000 K typically have M v ~ -2.0 (Mendez et al. 1992), which 
is more than 2 mag fainter than the ionizing stars of HII regions. 
Adopting a distance modulus of 25.6 mag for NGC 3 109 implies 
an apparent visual magnitude, V, brighter than about 21.3 mag 
for the ionizing stars in HII regions and fainter than about 23.3 
mag for the central stars of PNe. 

By measuring the visual magnitudes of the central stars of 
unresolved emitting objects, we can distinguish compact HII re- 
gions from PN candidates (see next section). Of course, it is pos- 
sible that a foreground or background star may appear projected 
on a PN candidate. In that case, our criterion would classify the 
object as an HII region rather than a PN, and spectroscopy is 
needed to classify the object properly. 

c. The r images 

The r filter has its maximum transmission at 6530 A and an 
FWHM of 810 A. It therefore includes the stellar flux and the 
nebular lines in this wavelength range. In particular, it includes 
Ha. Thus, we expect the [O in] -emitting objects to appear in 
the r images as well. Most of the objects do appear, in which 
cases they are most probably not background emission-line ob- 
jects. However, requiring detection in both [O m] and r discrim- 
inates against very faint (e.g., Magrini et al. 2001) but highly 
excited PNe that should have no stellar emission and for which 
the 5007/Ha ratio is large. Therefore, we will also consider with 
caution those candidates with faint [O m], but which do not ap- 
pear in the r image. 

d. The luminosity of recombination lines 

In addition to the previous criteria, it should be taken into ac- 
count that, often in HII regions, there is more than one ionizing 
star and these stars produce a larger number of ionizing photons, 
Qo, due to their larger radius. Typically log Qo > 48.0 (s~') for 
HII regions (Vacca et al. 1996) versus log Qo ~ 47.5 for PNe (us- 
ing evolutionary tracks by Blocker 1995). Optically thick com- 
pact HII regions should then be brighter than optically thick PNe 
in the hydrogen recombination lines. Furthermore, many PNe in 
the Milky Way are known to be optically thin, further weakening 
these lines. The Ha or H/3 fluxes (for instance) can be used for 
this criterion. 

e. The ionization degree or excitation class 

It is expected that, statistically, PNe will be more highly ex- 
cited than HII regions as a result of containing hotter ionizing 
stars. If a nebula shows He n 4686 emission whose intensity ex- 
ceeds a few percent of H/3, it is unquestionably a PN. Also, in- 
tensity ratios [O m] 5007 / H/3 larger than about 4 are found only 
in PNe. However, this criterion discriminates against PNe whose 
central stars have low T e jf (usually young PNe). Low-excitation 
PNe cannot be distinguished from HII regions on the basis of 
only their spectroscopic line ratios. 

From the above criteria, a and b are clearly the least ambigu- 
ous in distinguishing between compact HII regions and PN can- 
didates, when only imaging is available. The outcomes of our 
imaging analysis, upheld with the [O in] fluxes from our spec- 
troscopy, are presented in the following sections. 

3. Nebular and stellar magnitudes 

For all of our [O in] emitting objects, we have determined sev- 
eral instrumental magnitudes: the nebular magnitude m;(5007) 
from the continuum-subtracted on-band image, the stellar mag- 
nitude mi(star) from the off-band image, and m;(r) from the r 



image. Magnitudes were determined using the 'apphot' pack- 
age in IPvAF. Aperture photometry was performed by integrating 
the flux within a 5 pixel radius (equivalent to 1 arcsec) centered 
on the stellar position. The background was determined from an 
annulus 10-13 pixels from the star and subtracted from the flux. 
Special care was taken in crowded fields. Instrumental magni- 
tudes for each object are listed in Tables 1 and 2 for the PN 
candidates and the HII regions, respectively. 

3. 1 . Nebular [O m] 5007 magnitudes 

The nebular magnitude, m;(5007), was measured from the 
continuum- subtracted [O in] image, thereby excluding the emis- 
sion from the central star. The aperture employed includes all 
the flux for the nebulae with radii smaller than 1 arcsec, but it 
contains only a fraction of the flux for very extended objects. 
We find that the measured magnitudes for unresolved objects 
are well-correlated with the [O in] 5007 fluxes measured from 
our spectroscopy through the relation: 

log F(5007) = -0.40 mi(5007)- 8.32 . (1) 
This relation (shown in Fig. 1) has a small dispersion (<0.05 
mag) and a large correlation coefficient R 2 =0.997. This implies 
that the FORS 1 pre-imaging, reduced through the ESO pipeline 
(§2), is very consistent with the calibrated spectrophotometry. 
Thus, we are confident that the differences in instrumental mag- 
nitudes do represent flux ratios. 

3.2. Stellar magnitudes 

Apparent visual magnitudes for the central stars, m;(star), were 
measured in the off-band frame. A rough photometric calibra- 
tion using the same procedure as employed by Richer & McCall 
(1992) was used to relate these instrumental magnitudes to V 
magnitudes. That is, some stars from the list of Demers et al. 
(1985; hereafter DKI) with measured B and V magnitudes were 
used as local standards. Certainly, a calibration based upon ob- 
servations of photometric standards would have been more ap- 
propriate for a photometric calibration; but, as said in §2, our 
images were acquired as 'pre-imaging' for spectroscopy, so no 
standards were observed. 

Our off-band instrumental magnitudes follow a linear rela- 
tion with the V magnitudes from DKI, V(DKI), given by the 
expression: 

V(DKI)=0.99 mi(star) + 4.38 . (2) 
This relation is plotted in Fig.l. It has a dispersion <0.15 mag 
and a correlation coefficient R 2 =0.989. The good correlation of 
the two sets of magnitudes again indicates that differences in 
mj(star) indeed represent flux ratios, which allows us to assign a 
visual magnitude, V(DKI), to the central stars of our HII regions 
and candidate PNe. 

We also measured the instrumental r magnitude, m;(r). As 
already noted, m;(r) includes the stellar flux and the nebular lines 
in the r filter, particularly Ha. Consequently, objects with bright 
mi(r) magnitudes have very intense Ha or a bright star. 

3.3. PNe and HII regions in NGC 3109 

We have classified the diffuse and extended emission nebulae as 
HII regions. In the case of the compact objects, we separate HII 
regions from PN candidates using the magnitudes of their ioniz- 
ing stars. Considering that the visual magnitudes of the ionizing 
stars in HII regions should be brighter than about 21.4 mag (§2, 
b), equivalent to 17.2 mag in our system, we consider all objects 



4 



M. Pefia et al.: A catalog of planetary nebula candidates and HII regions in NGC 3109 



-13.0 



-13.5 



„ -14.0 
rv 
o 
o 

£ -14.5 



-15.0 



-15.5 



-16.0 



13.0 



14.0 



15.0 



16.0 
m r (5007) 



17.0 



18.0 



19.0 



19.0 



18.0 - 



q 17.0 



16.0 



15.0 



11.0 



11.5 



12.0 



12.5 



13.0 
mi(star) 



13.5 



14.0 



14.5 



15.0 



Fig. 1. The relations between our instrumental m;(5007) and the 
spectroscopic [O m] 5007 flux, F(5007), (upper panel) and be- 
tween nij(star) and V(DKI) for some stars (lower panel). 



with mj(star) brighter than 18 mag as HII regions. This limit- 
ing magnitude allows for an uncertainty of 0.3 mag in individual 
instrumental magnitudes (twice the dispersion of Eq. (2)) and 
an extinction of up to 0.5 mag, a value that amply includes the 
small color excess of E(B-V)= 0.04 mag (equivalent to Ay = 
0.13 mag, see §4) reported by Burstein & Heiles (1984) in the 
direction of NGC 3109 and any plausible internal extinction in 
the HII regions. We note that some ionizing stars, in both HII 
regions and PNe, could be binaries, in which case mj(star) and 
V(DKI) represent the fluxes of both stars. Considering the pos- 
sibility of binarity, our limiting magnitude for separating PNe 
from HII regions is then quite conservative. 

In Fig. 2, we present the color index m;(r)-m;(star) plotted 
as a function of mj(star) for HII regions (objects with mj(star) < 
18) and PN candidates (with m;(star) > 18). The colors of HII re- 
gions span the interval from -1.5 to 0.8 mag. For PN candidates 
the colors are bluer (more negative) than -1.3. The two types 
of nebulae are separated very nicely in this figure. This separa- 
tion occurs because the color index mi(r)-mi(star) corresponds 
roughly to the flux ratios (F(Har) + F r (star))/F,(star), which will 
be large for PNe while, at the same time, their central stars are 
faint. 

In total, we detected 20 PN candidates. Their characteris- 
tics are listed in Table 1, where we include their coordinates 
(columns 2 and 3) and magnitudes m;(5007), m;(star), V(DKI), 
and nii(r) (columns 4 to 7, respectively). Column 8 presents 
log F(5007), in erg crrT 2 s -1 , from our subsequent spectroscopy 
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Fig. 2. The color index mi(r)-mi(star) plotted against mi(star). 
Open squares represent HII regions (objects with m;(sfar) < 18 
mag). Black dots represent PN candidates. Since r includes the 
flux from the star and Ha, the color index is proportional to the 
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for stars with fainter m,.. 



(Pefia et al., in preparation), while column 9 lists the value of 
log F(5007) calculated from Eq. (1). In column 10, any previ- 
ous identifications are listed. As already reported in Pefia et al. 
(2006), we confirm as PNe only two of the seven PN candidates 
found by Richer & McCall (1992). 

Table 2 lists the characteristics of the HII regions. We present 
coordinates, diameters, magnitudes, and fluxes for 40 nebulae 
found in our continuum-subtracted [O m] image. This is not an 
exhaustive list of HII regions in NGC 3109, but it includes a 
large fraction of the [O m] 5007 emitting regions. The very low 
excitation nebulae without 5007 emission are obviously not in 
this list, and some detected but very faint HII regions and fila- 
ments were not included. The m;(star) and V(DKI) magnitudes 
measured for the central star of the HII regions (when there is 
an obvious central star) are presented in columns 6 and 7, re- 
spectively. In column 8 we list mi(r). The [O m] 5007 fluxes, 
F(5007), from spectroscopy and calculated via Eq. (1) are listed 
in columns 9 and 10, respectively. The other columns contain 
a description of the morphology (diffuse, compact, extended, 
complex, unresolved), a rough estimate of the stellar spectral 
type (see below), and previous identifications of the HII regions 
(columns 11-13, respectively). 

In Fig. 3 we present our [O ra] 5007 image of NGC 3109. 
In this figure, we show the distribution of PNe from Table 1 
and compact HII regions from Table 2. Our objects are spread 
throughout the galaxy, and it is evident that HII regions are rather 
concentrated towards the edge-on disk of the galaxy, while PNe 
are seen also above and below this structure, consistent with their 
belonging to an older stellar population. 

3.4. Individual PNe 

a. The brightest PN 

According to our selection criteria, our PN7 and PN 10 are 
the only bona fide PN candidates from the list reported by 
Richer & McCall (corresponding to their PN 1 and PN 6, re- 
spectively). PN 10 is a bright PN with no stellar continuum, but 
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Name 


R.A.* 


Dec* 


nii(5007) 


nii(star) c 


V(DKI) 




log rf 
F(5007) 


log' 
F(5007) 


PN 1 


10:02:41.57 


-26:08:40.7 


19.0 


>22 


>25 


> 22.5 




-15.91 


PN 2 


10:02:44.70 


-26:11:32.4 


19.7: 


21.0: 


25.1: 


17.8 




-16.21 


PN 3 


10:02:49.15 


-26:09:10.1 


16.4 


> 22 


>25 


18.8 


-14.88 


-14.89 


PN 4 


10:02:50.06 


-26:08:10.0 


17.3 


> 22 


>25 


> 22.5 


-15.25 


-15.23 


PN 5 


10:02:51.60 


-26:10:19.7 


18.7 


21.0: 


25.1: 


17.1 




-15.81 


PN 6 


10:02:57.20 


-26:08:55.8 


18.7 


19.5 


23.6 


17.9 




-15.81 


PN 7 


10:02:58.20 


-26:08:45.8 


14.4 


18.2 


22.4 


15.8 


-14.09 


-14.09 


PN 8 


10:02:59.52 


-26:08:13.9 


18.9 


20.0 


24.1 


17.4 




-15.89 


PN 9 


10:03:02.66 


-26:08:50.3 


16.4 


>22 


>25 


19.7 




-14.88 


PN 10 


10:03:03.91 


-26:09:42.3 


15.4 


> 22 


>25 


18.4 




-14.48 


PN 11 


10:03:05.04 


-26:10:54.2 


15.6 


> 22 


>25 


18.0 


-14.52 


-14.57 


PN 12 


10:03:07.99 


-26:07:57.8 


18.8 


19.40 


23.5 


16.5 




-15.21 


PN 13 


10:03:08.52 


-26:09:07.7 


17.2 


> 22 


> 25 


> 22.5 




-15.15 


PN 14 


10:03:10.34 


-26:08:06.1 


18.5 


19.5 


23.2 


17.1 


-15.78 


-15.71 


PN 15 


10:03:13.32 


-26:10:20.7 


19.3 


> 22 


>25 


> 22.5 




-16.04 


PN 16 


10:03:18.67 


-26:09:58.8 


17.0 


>22 


>25 


18.6 




-15.13 


PN 17 


10:03:18.74 


-26:10:05.4 


16.7 


18.8 


23.0 


17.5 


-14.96 


-15.01 


PN 18 


10:03:20.74 


-26:09:00.2 


17.5 


>22 


>25 


>22.5 




-15.32 


PN 19 


10:03:21.22 


-26:11:20.9 


19.1 


> 22 


> 25 


> 22.5 




-15.98 


PN20 


10:03:25.51 


-26:09:07.2 


17.0 


18.6 


22.8 


16.8 




-15.13 



PN1 RMc 
PN6 RMc 



(a) All these objects are unresolved. 

(b) Coordinates for equinox 2000. 

(c) Objects with m ; > 22 are not detected in the off-band image. Objects with m ; (r) > 22.5 are not detected in the r image. 

(d) Flux from spectroscopy, in erg cnT 2 s _1 . 

(e) Flux calculated from Eq. (1), in erg cnr 2 s' 1 . 

(f) RMc: Richer & McCall 1992. 




Fig. 3. Our [O m] on-band image, showing the whole galaxy. North is up, east is left. The full size of the image is 16.7 x 5.3 arcmin. 
PN candidates from Table 1 are marked with a circle. Compact HII regions from Table 2 are marked with a square (squares with an 
inner circle indicate the objects reported by Richer & McCall 1992). The HII regions are rather concentrated towards the edge-on 
disk of the galaxy, while PNe spread over a wider area. The large circles, from left to right, show the shells #1, 3, 5, and 9 studied 
by Hunter etal. (1993). 



PN7 has the brightest central star among our PN candidates, 
with V(DKT)=22.4 mag, i.e., brighter than expected for a sin- 
gle star (§2, criterion b). This nebula is also the brightest in [O 
m] 5007, one magnitude brighter than the next brightest object. 
We included it in the PN group due to its spectroscopy, which 



shows a very high 5007/H/? ratio (Pena et al., in preparation). In 
addition, it shows a [O m] 5007 line width, HWHM = 13 km 
s _1 (Richer 2006), which is significantly greater than the val- 
ues of HII regions. For comparison, the average [O m] 5007 line 
width (HWHM) of the five PN candidates from Richer & McCall 
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(1992) that we find are compact H II regions is 1 1 km s _1 (Richer 
2006; Richer et al. in preparation). Also in §4 we present further 
arguments for including it among the PN sample, 
b. The PNe with no or very faint r emission 

Six PN candidates in our sample were not detected in our r 
image. They are PN 1, PN4, PN 13, PN 15, PN 18, and PN 19. 
We determined an instrumental upper limit of 22.5 mag for their 
mi(r). Several of them are also very faint in [O m]. As mentioned 
in §2 (criterion c), these objects are considered PN candidates 
but some of them could be background sources. Most probably, 
though, they are high excitation PNe where [O m] 5007 is much 
stronger than Ha. From our spectroscopy, we have verified that 
at least PN 4 is a true PN (Pefia et al., in preparation). 

3.5. The stellar spectral types in HII regions 

On the V(DKI) magnitude scale, the central stars of HII regions 
show values in the range from 22.1 to 19.5, and a spectral type 
can be assigned to each them by adopting an absolute magnitude 
scale for massive stars (for instance Vacca et al. 1996 or Martins 
et al. 2005) and a distance modulus. We chose the distance mod- 
ulus, m-M= 25.62 mag (Minniti et al. 1999). This exercise was 
done only for those relatively small HII regions where there is 
an obvious central star. Extended nebulae (complex regions and 
loops) appear related to stellar clumps and it is very difficult to 
decide which ones are the ionizing stars; therefore for them we 
have not tried to identify the ionizing stars. Another problem is 
the assumption that these stars are single, which likely biases our 
classifications to earlier spectral types. In any case, these spec- 
tral types (column 9, Table 2) are very rough approximations 
that should be considered merely indicative. In addition, no red- 
dening has been assumed nor has any account been taken of any 
possible magnitude offset of the DKI calibration (Bresolin et al. 
1993 claimed a difference of about 0.3 - 0.4 mag). From this 
rough approximation we detect a few giants and supergiant stars 
(spectral type O III or O la), five main sequence stars between 
04 and 06, five 06-07, and 16 stars later than 07. 

4. The luminosity function for PNe 

From our [O m] fluxes, F(5007) (Eq. 1; column 9, Table 1), an 
apparent magnitude 1115007 can be obtained and the [O in] 5007 
luminosity function (PNLF) can be constructed. In Fig. 4, we 
present the cumulative luminosity function for our PN candi- 
dates (including PN7, see below). We chose this representa- 
tion because of the low number of PN candidates and because 
the peak magnitude of the luminosity function is poorly con- 
strained at the metallicity of NGC 3109 (Ciardullo et al. 2002). 
Therefore, we must fit the cumulative luminosity function using 
apparent magnitudes, instead of the usual absolute magnitudes. 
We convert the usual luminosity function 

N(M) oc exp(0.307M 5 oo7) (1 - exp(3(M 5 * 007 - M5007))) 
to a luminosity function in apparent magnitudes 

N(m) = N exp(-0.307//) exp(0.307m 5 (,o 7 ) (1 -exp(3(m; oo? - 
OT5007)) , 

where N is a normalization constant, /i is the apparent distance 
modulus, /i = 5 logd - 5 + A5007, M* 5im and m* QQ1 are the ab- 
solute and apparent peak magnitudes of the luminosity function, 
respectively, the latter defined via 

m 5 oo7 = -2.5 log F(5007) -13.74 
(Allen 1973; Jacoby 1989), and A 500 7 is the extinction at 5007 
A towards NGC 3109. For the extinction, we adopted the 
aforementioned color excess, E(B - V) = 0.04 mag, and the 



Table 3. Fit for the cumulative PNLF 



fit 




ra mo7 


A'exp(-0.307 y u) 


including 


PN7 


21.18 


2.62 x 10" J 


excluding 


PN7 


21.82 


2.77 x 10- 3 


LMC 




14.23 


0.287 


SMC 




14.82 


9.31 x 10~ 2 


M31 b 




20.20 


0.150 



(a) Based upon the sample of Jacoby et al. (1990) and 
considering the same limiting magnitude, 15.9 mag, for fitting. 

(b) Based upon the sample of Ciardullo et al. (2002) and 
fitting to the first three magnitudes of the luminosity function. 



Fitzpatrick (1999) reddening law, parametrized with a total-to- 
selective extinction of 3.041 . This parametrization delivers a true 
ratio of total-to-selective extinction of 3.07 when integrated over 
the spectrum of Vega (McCall 2004), which is the average value 
for the diffuse component of the interstellar medium of the Milky 
Way (McCall & Armour 2000). The extinction towards NGC 
3109 then amounts to A5007 = 0.13 mag. Our fit to the cumula- 
tive luminosity function returns the product N exp(-0.307//) and 
m* 5QQ1 as fitting parameters. Although this allows us to derive a 
normalization constant without knowledge of the peak magni- 
tude of the luminosity function, it obviously precludes using our 
peak magnitude to derive a distance to NGC 3109, as is usually 
done. 

Since the brightest PN (PN7 in Table 1, see §3.4) is 1 mag 
brighter than the next brightest object, we fit the cumulative lu- 
minosity function while both including and excluding this ob- 
ject. Both fits to the luminosity function are similarly good, so 
the luminosity function provides no grounds for excluding this 
object as a PN. The parameters for these fits are given in Table 
3, where we also provide similar fits to the cumulative lumi- 
nosity functions for PNe in the Magellanic Clouds and M31 
for comparison. Our apparent peak magnitudes for M3 1 and the 
Magellanic Clouds are within 0.12 mag of those found in the 
original studies (see Ciardullo et al. 1989 for M31), a difference 
that falls within the uncertainties in all of these studies, including 
ours. 

Considering the brightest PN candidate as a bona fide PN, 
we can use the true distance modulus of 25.62mag (Minniti et 
al. 1999) and the extinction A5oo7=0.13 to derive a peak absolute 
magnitude of the luminosity function of M* m)1 = -4.57 mag for 
NGC 3109. This value is surprisingly bright at the metallicity of 
this galaxy, compared with the value of -4.48 mag suggested by 
Ciardullo et al. (2002) from studies of several galaxies of higher 
metallicity. On the other hand, excluding the brightest PN, our 
fit to the luminosity function implies a peak absolute magnitude 
of only -3.93 mag, a value that is faint compared to the values of 
about -4.08 mag expected by Ciardullo et al. (2002) for galaxies 
with metallicities 12 -i-log O/H ~ 8.0. Obviously, one of the ben- 
efits of observing low metallicity galaxies like NGC 3109 will 
be a better understanding of the absolute luminosities of their 
planetary nebulae. 

Based upon the normalization factors in Table 3, NGC 3109 
has about the expected number of PNe. Adopting a true dis- 
tance modulus //()= 18.50 mag and A5007 = 0.36 mag for the LMC 
(Jacoby et al. 1990) and yUo=25.62mag and A5007= 0.13 mag for 
NGC 3109, we find N = 93.85 and N = 7.11 for the LMC and 
NGC 3109, respectively. If we adopt the absolute B magnitudes 
from Richer & McCall (1995), but correct them to the distance 
moduli already mentioned, the luminosity ratio is 11.69, very 
similar to the ratio of N values, 13.20. Therefore, NGC 3109 has 
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Fig. 4. Cumulative planetary nebula luminosity function from 
the data in Table 1 . The fit to the cumulative PNLF returns the 
product exp(-0.307ju) and m\ m (see text and Table 3). The 
LMC's luminosity function, scaled by a factor of 1/13.2 is also 
presented. 



about the same number of PNe per unit luminosity as the LMC. 
This is expected, since the colors of the LMC and NGC 3109 are 
similar (Buzzoni et al. 2006). 

In Fig. 4, we overlay the LMC's luminosity function, scaled 
by 1/13.2. Both fits (NGC 3109's and LMC's) show excellent 
agreement over the first 1.5-2.0 mag, the extent over which the 
LMC's luminosity function is believed to be complete (e.g., 
Richer et al. 1997). Since the normalization of the two fits to the 
luminosity function are similar (Table 2), regardless of whether 
PN7 is considered, the conclusion that NGC 3109 has about the 
number of PNe expected is independent of whether PN 7 may or 
may not be a true PN. 



5. Conclusions 

From [O m] 5007 on-band, off-band, and r imaging, we detected 
a large number of extended and compact emitting nebulae in 
NGC 3109. We have defined criteria to distinguish between PN 
candidates and compact HII regions, which allow us to identify 
20 PN candidates, for which we present coordinates, [O m] 5007 
magnitudes, and stellar magnitudes, when the central star is de- 
tected. We also present a catalog of 40 HII regions (extended 
and compact) that includes coordinates, diameters, [O m] 5007 
magnitudes, and instrumental magnitudes for the central stars. 
Apparent V magnitudes have been estimated from a rough pho- 
tometric calibration of the instrumental magnitudes. Based upon 
these V magnitudes, we assign a rough spectral class to the ion- 
izing stars of HII regions. 

The distribution of HII regions is concentrated towards the 
edge-on disk of the galaxy, while PN distribution is more ex- 
tended, covering a wider area, above and below the disk. This is 
consistent with the PNe belonging to an older stellar population. 

The cumulative PNLF for [O in] 5007 magnitudes was con- 
structed and a fit to this function was calculated. Comparing this 
fit with a fit to the available data for the LMC shows that NGC 



3109 has about the same number of PNe per unit luminosity as 
the LMC. By assuming a true distance modulus of 25.62 mag 
and an extinction A5oo7=0.13 mag for NGC 3109, a peak value 
of -4.57 mag is derived for the PNLF. This value is high com- 
pared both with the canonical -4.48 mag found for galaxies with 
higher metallicity and with the peak magnitude expected at the 
metallicity of NGC 3109. By studying more galaxies with metal- 
licities similar to NGC 3109's, a better understanding of the ab- 
solute luminosities of low metallicity PNe will be acquired. 
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Table 2. Positions and characteristics of HII regions in NGC 3109 
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HII 27 


1 A AO 1 O O A 

10:03:12.84 


O ^ A A C A A 

-26:09:54.9 


2.0 


16.8 


17.9 


22.1 


17.2 




-15.05 


din 


/~\ O r\A T 7 

08-09 V 




TTTT lO 

Hll 28 


1 A. AO .10 OO 

10:03:13.32 


O/C.AO.CA A 

-26:08:50.0 


O A 

LA 


1 C 

15.6 


1 H A 

17.0 


21.2 


1 *7 1 

17.1 




1/1 CO 

-14.58 


din 


t~\o \ J 

08 V 




HII 29 


1 A AO 1 ^7 /I O 

10:03:17.48 


i"\ /" 1 a A A C 

-26:10:00.5 


6.0 


16.7 


16.9 


21.1 


17.1 




-14.99 


j ■ re 

din 


07-08 V 


II1T1 A n»/i TT O 

HII 10 RMc, Ho2 


TTTT OA 

Hll 30 


1 A. AO . 1 1 7 1 

10:03:17.71 


O/C.AA.CA /I 

-26:09:59.4 


A O 

0.8 


1 £ A 

16. 


1 *7 C 

17.5 


O 1 H 

21.7 


17.6 




— 14.74 


comp 


08-09 V 


near Hll 10 KMc 


HII 31 


1 A AO OO AA 

10:03:22.00 


O tl AO CO A 

-26:08:53.4 


o c 

2.5 


1 c o 

15.8 


1 *7 tl 

17.6 


21.8 


1 "i ti 
17.6 




1/1 ^ c 

-14.65 


j ■ re 

din 


08-09 V 


D1H6 


HII 32 


1 A AO OO AO 

10:03:22.92 


O tl 1 A OO /I 

-26:10:2 7 .4 


1.5 


15.6 


16.5 


20.7 


16.2 


-14.58 


-14.55 


comp 


07 V 




HII 33 


1 A AO OC CA 

10:03:25.50 


i"\ /" 1 A C /I ~ 7 

-26:10:54.7 


18.0 


17.5 


16.3 


20.5 


16.9 




-15.33 


ext 


05-06 V 




T TTT 1 A 

Hll 34 


1A.AO.OC OA 

10:03:25.80 


O/C.AA.AO A 

-26:09:08.9 


A A 

0.9 


17.4 


1 H A 

17.0 


O 1 O 

21. 2 


i ti n 

16.7 


1 C OO 

-15.23 


1 C OA 

-15.29 


comp 


07 V 




HTT 

nil jj 


1 n-n^-?7 1 Q 

L\J.\JJ.Z> 1 . 17 


?fv 1 0-77 4 




1 7 7 


17 ^ 




17 S 




— 1 s 




f«lii m n 


Ha1 Hh1 
nu i , nui 


HII 36 


10:03:27.90 


-26:09:36.4 


2.5 


17.3 


17.4 


21.6 


16.4 




-15.25 


diff 


08-09 V 




HII 37 


10:03:27.96 


-26:09:19.9 


0.9 


13.6 


16.7 


20.9 


15.2 




-13.75 


comp 


07-08 V 




HII 38 


10:03:28.00 


-26:10:22.8 


3.0 


16.9 


16.8 


21.0 


17.6 




-15.09 


diff 


07-08 V 


F5H1 


HII 39 


10:03:28.08 


-26:09:14.8 


3.0 


17.5 


17.3 


21.5 


17.8 




-15.33 


diff 


08-09 V 


F5H2 


HII 40 


10:03:28.51 


-26:09:21.8 


0.8 


13.2 


15.4 


19.6 


14.4 


-13.55 


-13.57 


comp 


clump 


near F5H4,H3,H5 



(a) V stellar magnitude of the ionizing star, as derived relative to DKI photometry. 

(b) Flux in erg cm~ 2 s~', column 9: from spectroscopy; column 10: from expression (1). 

(c) Nomenclature: comp: compact, cmplx: complex, diff: diffuse, ext: extended, unres: unresolved. 

(d) References: RMc, Richer & McCall 1992; Ho: Hodge 1969; Hu: Hunter et al 1993; F3H1, F5H4, D1H6, ... : Bresolin et al. 1993. 



